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Abstract 

A methodology for simulation of molecular mixing, and the 
resulting velocity and temperature fields has been developed. 
The ideas are applied to the flow conditions present in the 
NASA Lewis Research Center Planar Reacting Shear Layer 
(PRSL) facility, and results compared to experimental data. A 
gaussian transverse turbulent velocity distribution is used in 
conjunction with a linearly increasing time scale to describe 
the mixing of different regions of the flow. Equilibrium 
reaction calculations are then performed on the mix to arrive 
at a new species composition and temperature. Velocities are 
determined through summation of momentum contributions. 
The analysis indicates a combustion efficiency of the order of 
80 percent for the reacting mixing layer, and a turbulent 
Schmidt number of 2/3 . The success of the model is attributed 
to the simulation of large-scale transport of fluid. The favor- 
able comparison shows that a relatively quick and simple PC 
calculation is capable of simulating the basic flow structure in 
the reacting and nonreacting shear layer present in the facility 
given basic assumptions about turbulence properties. 

Nomenclature 

b T local mixing layer thermal width, m 

b y local mixing layer velocity width, m 

db/dX mixing layer growth rate 
I integral time scale, s 

U streamwise velocity, m/s 

U mean streamwise velocity 

V' transverse velocity, m/s 

V' ms root mean square transverse turbulent velocity, m/s 
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y mean transverse velocity, m/s 
x streamwise distance, m 

y transverse distance, m 

p autocorrelation coefficient 

t time delay, s 


Introduction and Previous Work 

This report summarizes the theory behind mixing layer 
simulations created through the use of a mixing simulation 
program and the NASA Lewis CETPC Chemical Equilibrium 
Composition program 1 , and presents some initial results. The 
mixing program, entitled MYV.FOR, was created to simulate 
the velocity, thermal, and species distributions in the NASA 
Lewis Research Center Planar Reacting Shear Layer (PRSL) 
facility. 2 In evaluation of the modeling, experimental thermo- 
couple and Laser Doppler Velocimetry (LDV) data is used. 3,4 
MYV.FOR is not a fluid dynamics program, but rather a 
methodology for simulation of molecular mixing. A transport 
model based on transverse turbulent velocity profiles and 
integral time scales is used to simulate the mixing of mol- 
ecules in the flowfield. Through conservation of momentum, 
a new velocity distribution is determined. The mix is assumed 
to react and convect to a downstream location. The new 
composition is the initial condition for subsequent mixing, 
reaction, and convection at the downstream location. 

Motivation for the current effort lies in the large amount of 
time required to generate a full Navier-Stokes reacting flow 
simulation of the mixing layer. It is desired to develop and 
demonstrate a simple model which can run in minutes on a 
typical desktop computer. Advantages include both the time- 
liness of the results, and the insight into the turbulent transport 
gained. 

Previous work has been done for the reacting shear layer 
configuration, both experimentally and computationally, by 
many investigators. Because the emphasis of this work is on 
the simplicity and speed of the model, only simple, 


empirically based works are mentioned here. Konrad 5 mea- 
sured the probability density function (PDF) of molecular 
constituents in a mixing layer formed from two dissimilar 
streams. He computed a product PDF based on the two 
freestream reactant PDFs for a theoretical reacting mixing 
layer. For comparison purposes, he compared that result with 
that resulting from using only mean reactant profiles, which is 
similar to the current effort. He found that use of the PDF 
model lowered temperatures and product fraction signifi- 
cantly below the equilibrium calculation. 

Wallace 6 used Konrad’s nonreacting PDF’s in a similar 
calculation, with the goal of simulating his experimentally 
obtained reacting mixing layer thermal profiles. The compu- 
tational profiles were, in general, higher in temperature, 
shifted in position, and over-predicted diffusion. 

In contrast, the current approach uses the transverse veloc- 
ity PDF in an attempt to model time averaged momentum 
diffusion. An assumed Schmidt number is used to model 
molecular mixing. Then, an equilibrium calculation is per- 
formed to arrive at post-reaction molecular compositions, 
temperatures, and densities. Momentum mixing is performed 
as well, and mean streamwise velocity profiles generated. 

Computational Domain and Experimental Conditions 

A schematic of the experimental facility is shown in Fig. 1 . 
The upper stream typically flows at a Mach number of 0.4. For 
nonreacting tests, air is used in both streams, while for the 
reacting mixing layer studies, the upper stream consists of 
hydrogen diluted with nitrogen. The lower stream for both 
reacting and nonreacting studies is heated air flowing at a 
Mach number of 0.7. The test section is approximately 40 cm 
in the streamwise direction, 1 0 cm in the vertical direction, and 
20 cm in width. 


The facility was designed to investigate the compressible 
flow region of high subsonic Mach numbers. The convective 
Mach number, defined by the velocity differential between the 
convective velocity and one freestream, and the speed of 
sound in that freestream, is valued at approximately 0.28. The 
Reynolds number based on the velocity differential, the aver- 
age viscosity, and the layer width at 300 mm downstream of 
the knife edge, is approximately 1.8xl0 5 . For further discus- 
sion of facility design and flow conditions, as well as a detailed 
summary of experimental results, see Refs. 2 to 4. 

For computational purposes, the PRSL test section is di- 
vided into 400 transverse elements. Experimentally, the 
reacting mixing layer is found to impinge the top of the test 
section at far downstream positions, due to the large growth 
angle and upward slant of the layer. The computational 
method simulates transverse fluid movement from each 
freestream into the mixing layer. Thus, the computational 
domain extends beyond the dimensions of the actual test 
section to avoid difficulties with the statistical mixing proce- 
dure. Because of the large transverse influence distance at 
these downstream locations, the computational domain is 
actually four times the height of the test section. Each compu- 
tational element is initially 1 .01 6 mm in height. The lengths of 
the computational elements are determined by the streamwise 
convection distance which will be discussed further. 

Since the facility is designed to impose a neutral axial 
pressure gradient through adjustable walls, and experimental 
data shows that this is reasonably accomplished, constant 
pressure combustion calculations are used. Thus, density 
changes considerably with heat release, by more than a factor 
of three near stoichiometry. Because no axial acceleration of 
the fluid is observed, and the fluid is constrained to a constant 
width in the cross-stream direction, all growth necessary for 
the conservation of mass is assumed to take place in the 
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Figure 1 . — Planar reacting shear layer facility schematic. 
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transverse direction; that is, the initial height of the fluid 
elements described above is allowed to increase to account for 
fluid expansion. 

Computational Method 

The method can be summarized as a statistical molecular 
mixing simulation followed by equilibrium combustion and 
convection downstream of reaction products. The important 
parameters governing the mixing are the transverse turbulent 
velocity probability density function (PDF) and the integral 
time scale / of the turbulent flowfield. Each fluid element is 
allowed to disperse and travel in the transverse direction. The 
relative amount of each element travelling at a given trans- 
verse velocity is determined by the PDF. The transverse 
velocity, coupled with a time scale, determine the cross- 
stream distance travelled by all portions of each fluid element. 
The same time scale is used in conjunction with the average 
streamwise velocity to determine downstream convection 
distance. 

For the transverse turbulent velocity a gaussian PDF is used 
throughout the presented simulations. The selection of this 
PDF function is discussed in the Computational Parameters 
section. Thus, for zero mean transverse velocity during a 
mixing simulation, the majority of a given fluid element will 
have zero cross-stream velocity, while decreasing amounts 
will have both positive and negative transverse velocities of 
increasing magnitude. 

The integral time scale I should be thought of as a ‘coher- 
ence time* for fluid eddies during travel. Each PDF weighted 


portion of a fluid element is assumed to travel for time I before 
mixing with other arriving fluid elements. Empirically, we 
expect the integral time scale for the mixing layer to increase 
linearly with streamwise distance. Selection of /, and the effect 
on computational layer growth rate, is discussed in the Com- 
putational Parameters section. 

Computationally, we look at the contributions to a point of 
interest from all other elements in the flow. Thus, we shift here 
from a description of the dispersion and travel of a given fluid 
element to a calculation of contributions from all elements to 
a point of interest. Figure 2 shows a schematic of three 
contributing elements to a computational point and their travel 
distances. Elements A,B, and C are transverse distances y(A), 
y(B), and y(C) from the element of interest D. There exists a 
local integral time scale / over which fluid may travel. The 
combination of transverse distance and travel time produce a 
required transverse velocity for contribution of A, B, or C to 
D as 



where y' is the transverse velocity and the index represents 

transverse fluid elements. Using the gaussian cross-stream 
velocity PDF for A, B, and C, we can determine the portion of 
each element which will possess the required transverse 
velocity to travel to D in time I through 

1 s {-v’-vf/w£s 

V 27tV'rms 


Mixing volume 
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The product of the number density in molecules/cc and vol- 
ume of the contributing element are multiplied by this weight- 
ing to determine the number of molecules arriving at D. The 
enthalpy of each contributing fluid element is summed to 
arrive at a total number and enthalpy of each molecule at D. All 
species are then assumed to mix and react to equilibrium. The 
equilibrium composition and enthalpy is determined through 
the use of the NASA Lewis CETPC code. 1 The momentum 
contributions are similarly summed, and used to arrive at a 
mixture streamwise velocity. Also noted on Fig. 2 is the 
streamwise convection distance of the resulting equilibrium 
composition, given by 

x = UxI 

where the average of the two freestream velocities is used. 

Inherent in the methodology are many simplifications, 
necessary to provide the code with if s desired simplicity and 
speed. One simplification is the time-averaged nature of the 
calculation of contributions. This description assumes that all 
elements at the streamwise calculation location contribute to 
the element of interest at all times. Thus, the time-dependent 
nature of turbulent mixing is not modeled. It is also assumed 
that upon reaching the central element, all contributions mix 
well to equilibrium. To provide a correction for these two 
assumptions, an efficiency term is introduced. This idea of a 
time-averaged efficiency was developed while observing Laser 
Induced Fluorescence (LIF) images of the reacting shear 
layer. It was noted that combustion products were present in 
small pockets, or eddies, sporadically throughout the mixing 
layer. By extending the physical observations to a time aver- 
aging of post-combustion temperature and freestream tem- 
peratures, the experimentally observed thermal profiles for 
the reacting shear layer were very closely reproduced. 

This efficiency term is used as the fraction of time that the 
calculated molecular distributions are indeed present, and 
react to equilibrium. For the remaining time, each of the 
contributing elements are assumed to be present according to 
their statistically weighted contribution. The calculated tem- 
perature, then, is a time average of the reacted mix and the 
component elements. The calculated width of the reacted 
element is corrected by the ratio of the averaged and post- 
reaction temperatures. Following convection, subsequent 
mixing calculations are based on the equilibrium composition. 
In the results presented herein, an efficiency of 80 percent has 
been used. This value produces an approximate match be- 
tween the peak experimental and calculated temperatures. 

Due to the nature of the computational method, profiles are 
generated only at certain streamwise positions as determined 
by the product of average streamwise velocity and local 
integral time scale I. Experimental data, too, exists only at 
specific streamwise positions. To allow comparisons between 
experimental and computational results, an interpolation rou- 


tine for the calculated flowfield was developed to provide 
profiles at desired streamwise positions. For the nonreacting 
case, this consisted of nondimensionalizing the transverse 
velocity and thermal profiles into two curves, both similar to 
the error function. These nondimensionalized profiles were 
found to be self-similar at all streamwise locations, allowing 
interpolation at intermediate stations. Linear fits to transverse 
width and centerline position were used to determine required 
fitting parameters at any desired streamwise position, and 
profiles generated. For the reacting case, the additional com- 
plication of heat release is present. It was found that by 
modeling a nonreacting thermal layer afterthe nondimensional 
velocity layer, fitting a nonreacting thermal profile, and 
subtracting it from the computed reacting thermal profiles, a 
somewhat bell-shaped curve of temperature rise was formed. 
Dividing this curve by it ’ s maximum value led to a self-similar 
family of nondimensional temperature rise curves for all 
streamwise positions. By fitting lines to the maximum tem- 
perature rise, it’ s transverse position, and the temperature rise 
width as a function of streamwise distance, reacting thermal 
profiles were generated at desired streamwise locations. 

Throughout model development, a measure of the growth 
rate of the calculated mixing layers was needed for compari- 
son with experimental data. The 10 percent delta thicknesses 
bj and by were used for thermal and velocity thicknesses 
respectively. The upper and lower bounds of the mixing layer 
are denoted as the point where the local condition departs from 
the near freestream condition by 10 percent of the freestream 
velocity or temperature differential. The distance between 
these two bounds is the local width b of the mixing layer. The 
growth rate is the width b divided by the streamwise distance 
X from the splitter plate tip, as 

db__b_ 

dx~ X 

Computational Parameters 

The functional form for the integral time scale / was chosen 
through comparison of experimental growth rates with com- 
putational results. It was determined that a constant value of I 
produced a narrow computational mixing layer. To obtain 
growth rates similar to those observed experimentally, the 
time scale had to increase with downstream distance. A 
linearly increasing time scale is in accordance with intuitive 
expectations. The effect of the slope of I on mixing layer 
growth rate was investigated for a linear formulation. For the 
nonreacting case, the mixing layer velocity width b v was 
found to vary as the slope of I raised to the 1/4 power. The 
effect of a large nonzero intercept for the I versus X curve is 
to provide a rapid growth region at the beginning of the mixing 
layer. The subsequent growth rate is not affected. For compu- 
tational purposes, a nonzero intercept is required. The inter- 
cept used, 0.01 ms, is small enough to provide an effectively 
zero initial width of the layer. 
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Experimental determination of the flowfield integral time 
scales yields varying values in both the streamwise and 
transverse directions. The integral time scale is the integral of 
the autocorrelation coefficient, as 


/ = J p{T)dx 
0 

where 

p(T) = =L= x [u(t) x u(t + t)] 
u 2 (t) 

In general, values of 7 peaked within the mixing layer region 
and increased with downstream distance. The current resolu- 
tion of experimental data, however, is insufficient to define the 
form of 7 in the transverse direction accurately. To keep the 
computational model simple, it was desired to use a constant 
value in the cross-stream direction. Maximum experimental 
values of 7 divided by the position X of measurement yield a 
slope of the 7 versus X curve of the order of 0.006 s/m; this is 
the value used throughout all simulations. Use of the smaller 
7 values present in the freestream results in a narrow compu- 
tational layer. Insertion of a transverse 7 profile in keeping 
with experimental observation would have introduced unde- 
sirable complications into the model. 

As mentioned in the Computational Model section, a 
gaussian PDF was used to describe the transverse turbulence 
profile. Analysis of experimental data 3,4 indicates this is a 
reasonable approximation. As in the case of the integral time 
scale, actual experimental data for the transverse turbulence 


PDF varies in both the streamwise and cross-stream direc- 
tions. Again, undesirable complications arise from the use of 
nonconstant values. The relationship between the rms 
value of transverse velocity and the mixing layer growth rate 
was investigated. The velocity layer growth rate of the 
nonreacting mixing layer was found to be given approxi- 
mately by 


for the formulation of 7 di scussed previously . A corresponding 
simple expression for reacting mixing layers is not available 
due to the expansion that occurs in accordance with heat 
release. The factor of three in the above expression arises from 
the gaussian transverse turbulent profile. Three times the root 
mean square value encompasses 99.7 percent of the area under 
a gaussian curve, and very little fluid with velocity magnitudes 
in excess of this figure move throughout the flowfield in this 
formulation. 

By equating the above expression with the experimental 
nonreacting mixing layer growth rate, a value of V' ms - 
14.6 m/s was found to reproduce the nonreacting velocity 
layer growthrate.Experimental values of V' ms at all streamwise 
measurement locations nondimensionalized by average 
streamwise velocity are shown in Fig. 3 along with the 
empirically determined computational value. This value is on 
the high end of experimentally observed transverse velocity 
rms values. Use of an average value results in a computational 
mixing layer more narrow than the experimentally observed 
result. Thus, as in the case of the integral time scale, the model 
requires the use of the maximum value present in the flowfield 
to reproduce the experimental result. 
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Table 1. - Flow Conditions And Statistical Data For Simulations 


Input Values 

Reacting flow 

Nonreacting flow 

Lower duct air flow, kg/s 

1.94 

1.87 

Upper duct N 2 flow, kg/s 

1.00 

1.62 

Upper duct H 2 flow, kg/s 

0.032 

0.000 

Lower duct velocity, m/s 

375 

380 

Upper duct velocity, m/s 

135 

134 

Lower duct static temperature, K 

740 

758 

Upper duct static temperature, K 

322 

335 

Static pressure, Pa 

106 300 

100 000 

Average transverse velocity, m/s 

12 

0 

RMS transverse velocity, m/s 

14.6 

14.6 

Intergral time scale, s (X in m) 

0.0001 + 0.006*X 

0.0001 + 0.006*X 


A value for the Prandtl number, or ratio of momentum 
transport to thermal transport, was chosen based on the rela- 
tive widths of computed nonreacting thermal and velocity 
mixing layers as 2/3. The same value was assumed for the 
Schmidt number. This value was implemented as a multiplier 
of the rms transverse velocity value to effect a change in 
molecular transport. 

Table 1 summarizes flow conditions for the two simula- 
tions. In the nonreacting case, pure nitrogen was used to 
simulate the air stream. The Results section presents compari- 
sons between experimental and computational profiles. No- 
ticeable in many of the figures are discrepancies between 
computational and experimental freestream flow variables. 


These arise from the limitations of both experiment and 
computation. Computationally, we have constant freestream 
temperatures and velocities. Experimentally, freestream tem- 
peratures and velocities tend to drop somewhat with streamwise 
position. Thus, it is impossible to match freestream conditions 
at all positions of comparison. Freestream velocities were 
chosen to nearly match the majority of the data, along with 
temperatures. Due to this small departure from specified 
experimental inlet conditions 3,4 , adjustments to the freestream 
static pressure of the order of one percent were necessary to 
match the experimentally specified inlet mass flow rates. 

The values of V' ms , /, and turbulent Schmidt number 
determined as previously discussed for the nonreacting mix- 
ing layer were applied to both flowfields. These are the most 
significant variables in determining mixing layer growth rate. 
The large differences in growth rate observed computationally 
and presented in the results section are due to the expansion 
which occurs under the constant pressure combustion for the 
reacting mixing layer. 

Results 

Figure 4 shows the calculated mixing layer temperature and 
velocity contours for the reacting shear layer simulation. Due 
to the positive mean transverse velocity, the mixing layer 
moves upward in the duct. The mixing layer widths increase 
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Figure 4. — Computed reacting mixing layer temperature and velocity contours. 


6 




4 


linearly with streamwise distance. Similar flowfield visual- 
izations for the nonreacting case show the centerline of the 
mixing layer to be below the duct midline, in spite of the zero 
mean transverse velocity. This is due to the momentum and 
enthalpic mixing of the low density high speed lower stream, 
and high density low speed upper stream. This observation 
agrees with the experimentally observed result. As will be 
seen in data to be presented, the reacting layer is found to be 
much wider than the nonreacting layer due to the expansion 
of the constant-pressure combustion products. This, too, is in 
qualitative agreement with experimental data. 

Shown in Fig, 5 is a comparison between experimental and 
calculated temperatures at a position 150 mm downstream 
from the splitter plate tip. For the nonreacting case, the 
calculation reproduces the growth rate of the experimentally 
measured thermal mixing layer well. As this case was utilized 
to determine the Prandtl and Schmidt numbers to use for the 
simulation, the agreement in shear layer width is expected. 
The computational profile is centered approximately 2 mm 
high relative to the experimental data, however. The Schmidt 
number determined through the nonreacting thermal width 
was used in the calculation of the reacting thermal profile, also 
shown in Fig. 5. The experimental reacting profile is similar 
in shape to the experimental result. However, the experiment 
profile shows a thicker region of high temperature gas. As in 
the nonreacting case, the vertical position of the experimental 
thermal layer is poorly predicted. The prediction of vertical 
position of the thermal and velocity mixing layers will be 


discussed following presentation of all comparisons between 
experimental and computational data. 

Figure 6 shows the temperature profile data 300 mm down- 
stream from the splitter plate tip. As in the case of the 1 50 mm 
data, good agreement is seen between the computational and 
experimental results in width of the nonreacting thermal 
mixing layer. Again, the vertical positioning is not repro- 
duced. For the reacting shear layer, the comparison between 
the computational and experimental result is excellent. The 
general shape of the thermal layer, as well as the width of the 
high temperature region, is well-matched. 

Figure 7 shows the calculated and experimentally obtained 
velocity profiles for nonreacting flow. While the vertical 
positioning of the velocity mixing layer is again poorly mod- 
eled, the diffusion is closely matched for both the X = 1 50 mm 
and X = 300 mm streamwise positions. As in the case of the 
nonreacting temperature profile width, this agreement is ex- 
pected due to the empirical nature of the formulation for the 
nonreacting case. The nonreacting mixing layer widths were 
used in conjunction with the growth rate formula presented 
previously to arrive at the value of the transverse turbulent root 
mean square velocity used in computations. 

Applying the values obtained through the nonreacting 
studies, the model produces the reacting layer velocity pro- 
files shown in Fig. 8. Here, we see excellent agreement 
between experimentally measured velocity profiles and the 
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computational result for shear layer width at both positions of 
comparison. The X = 150 mm streamwise position shows a 
slight discrepancy in the vertical placement of the mixing 
layer between the computational and experimental result. The 
X = 300 mm data, however, is well-modeled both in layer 
width and vertical position. 


Discussion 

As mentioned previously in presentation of the computa- 
tional results, the good agreement between computed and 
measured thermal and velocity mixing layer widths is ex- 
pected for the nonreacting simulations. This is due to the 
empirical method used to arrive at a formulation of /, the value 
of the rms transverse velocity, and the Schmidt number. These 
three quantities are the significant flowfield variables in 
control of mixing layer width in this model. The test of the 
method lies in whether the shape of the nonreacting profiles 
are matched, and whether the flowfield parameters obtained 
through analysis of nonreacting mixing layers can success- 
fully model the reacting shear layer. 

For evaluation of the computational work with the experi- 
mental data, there are three comparisons available. One is the 
general shape of the contours. The experimental nonreacting 
temperature data, and velocimetry data from both reacting and 
nonreacting cases, has been shown 3,4 to collapse very nicely 


onto an error function curve. The calculated profiles differ 
slightly from the error function due to the momentum and 
enthalpic mixing. Collapsed thermal and velocity nonreacting 
profiles are shown in Fig. 9, along with the error function. If 
simple averages are used of temperatures or velocities, the 
gaussian mixing used will generate error function profiles as 
a result. However, the weighting afforded by streams of 
differing densities and enthalpies results in computational 
profiles that differ somewhat from that of the error function. 
This points out a limitation of the computational method; the 
shapes of the profiles generated are fundamentally somewhat 
different than those observed experimentally. This limitation 
is observable in the computational profiles presented previ- 
ously as a more gentle comer on the upper side of the velocity 
profiles and nonreacting thermal profiles as compared to the 
experimental data. The discrepancy, however, is not large. 

A second criterion is the modeling of diffusion rate for the 
shear layer with heat release. Looking at the comparisons 
between widths of experimental and computational velocity 
layers, we see that the growth rates are very well modeled with 
the assumed integral time scale and transverse velocity PDF 
for the reacting mixing layer. As mentioned previously, the 
large relative velocity width b v of the reacting mixing layer 
relative to the nonreacting case is due entirely to thermal 
expansion in the reacting case. All flow variables which 
contribute significantly to growth rate were held constant 
between the two cases. The thermal layer at 300 mm 
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Figure 9.— Computational non-dimensionalized profiles and the 
error function. 

downstream of the splitter plate is well modeled in shape and 
width. From these comparisons, it appears that the diffusion, 
heat release, and subsequent layer expansion is well modeled. 

The discrepancy between computational and experimental 
150 mm thermal data is difficult to explain with this in mind. 
The shape of the reacting thermal layer here seems to be the 
problem as much as the width. Experimentally, the mixing 
layer region appears to have expanded to a larger degree than 
in the computational result. An investigation of the effect of 
Schmidt number and efficiency on the reacting thermal layer 
shape provided no means of matching the observed profile. 
Simulation of a higher diffusion rate, through use of a smaller 
Schmidt number, tends to increase the width of the wings of 
the calculated layer more than the middle. Increasing the 
efficiency increases heat release, and so results in greater 
expansion; however, the resulting temperature increase at all 
points in the flowfield results in a computational profile that 
simply overshoots the experimental data everywhere. Thus, at 
this time, the discrepancy is unexplained. Further examination 
of additional experimental thermocouple data may shed some 
light on the problem. 

A third criterion is the vertical placement of the mixing 
layer, which is rather poorly modeled in many cases. This is to 
be expected to a certain degree. The code is incapable of 
simulating any fluid dynamic effects which might result in 
curvature of the mixing layer, or movement in any other than 


a constant angle with respect to test section centerline. How- 
ever, the mismatch may be partially attributed to limitations of 
the experimental data. The experimental thermal profiles 
presented for a given condition were obtained on different 
dates, and none in conjunction with the velocimetry data 
presented. Thus, the vertical placement, controlled through 
the mean transverse velocity, is effectively a guess for the 
thermal profiles. Data has shown that the mixing layer posi- 
tion in the vertical direction changes considerably from test to 
test. The phenomenon behind this placement of the layer is not 
yet understood. The current feeling is that nonzero mean 
vertical velocities are the cause; but, the source of the nonzero 
mean velocities is not yet understood. The vitiating torch 3 used 
to ignite the mixing layer is one possible source of the vertical 
velocity component. Inconsistencies between mixing layer 
position in the transverse direction at the two streamwise 
positions are therefore not regarded as a serious limitation of 
the computational method at this time. Further experimental 
investigations will continue to define the computational bound- 
ary conditions. 

The discrepancies in the position of the velocity mixing 
layer, too, are not regarded as a serious deficiency. It has been 
determined that the test section moves several millimeters 
both vertically and horizontally during testing due to thermal 
expansion. Feedback sensors to account for this expansion 
have recently been installed on the test rig. Additional data 
with the origin corrected should remove any ambiguity from 
the vertical and horizontal distance between the splitter plate 
tip and the measurement position, allowing continued code 
evaluation and refinement. 

The success of the code in simulating both the nonreacting 
and reacting thermal layers through simulation of molecular 
transport is due to the fact that large scale movements of gases 
are allowed in the model. The large influence distance dis- 
cussed previously allows reactants to penetrate deeply into the 
computational mixing layer, resulting in a large heat release. 
The success of the model indicates that the simulation of large- 
scale transport of molecules is necessary to reproduce this 
highly turbulent flowfield. Molecular diffusion, as opposed to 
turbulent mixing, is not modeled here except in the form of the 
inefficiency discussed previously. Increasing the efficiency, 
and hence gaseous diffusion rates, does not reproduce the 
wide thermal layers without the mixing-based handling of the 
Schmidt number. 


Code Application 

The results presented herein are the only ones generated to 
date by this code. The keep the code as general and simple as 
possible, the values root mean square transverse turbulent 
velocity as well as integral time scale were held constant 
throughout the flowfield. With these simplifications, it was 
found that the use of the maximum values of turbulence 
parameters present in the shear layer was necessary to repro- 
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duce the experimental thermal and velocity profiles. A desir- 
able extension of the present work is to apply the code to 
alternate geometries and flow conditions, and evaluate the 
performance. 

Suitable for testing are high Reynolds number annular or 
planar mixing layers. Recommended flowfield variables, based 
on the current result, are a root mean square transverse 
velocity equal to 5.7 percent of the mean freestream velocity; 
a turbulent Schmidt number of 2/3; and an integral time scale 
formulation of / = 0.00001+.006*X in units of seconds, with 
streamwise distance X in meters. A mixing and combustion 
efficiency of 80 percent could be used initially, with correc- 
tion made to match peak experimental temperatures if such 
data were available. 


Conclusions 

A simple methodology has been developed for the simula- 
tion of the PRSL facility reacting flowfields. Mixing rates 
between the two streams are simulated relatively accurately, 
with the momentum mixing fitting experimental data more 
closely than the thermal. Mixing layer centerline position and 
curvature are not well modeled, but the discrepancies are not 
viewed as serious limitations of the model at this time due to 
experimental uncertainties. The simulations suggest a mixing 
layer efficiency of approximately 80 percent and a turbulent 
Schmidt number of 2/3 . Identical turbulent flowfield variables 
for both cases successfully reproduce experimental data, 
suggesting that the fundamental properties of the flowfields 
are similar. The success of the model indicates that simulation 
of large-scale transport is necessary to reproduce the observed 
properties of this highly turbulent mixing layer. The calcula- 
tion is simple enough to run in a few minutes on a desktop 
computer, and shows that a simple model based on bulk 
transverse turbulent transport can effectively reproduce ther- 
mal and velocity experimental data for this facility. The 
formulation has not yet been applied to any other facility or 
flow geometries, and should be compared with additional 
experimental data before it is applied to an untested flowfield 
to generate profiles. 
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